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We present the doping dependence of the spectral functions, density of states and low frequency 
behavior of the self-energy for the 2D Hubbard model on the basis of our recently developed theory for 
the Hubbard model. Strong 2D critical spin fluctuations dominating near half-filling are completely 
destroyed by 13 and 20 % doping concentrations for U = 4 and 8, respectively. Below these 
concentrations the imaginary part of the self-energy vanishes quadratically in frequency near the 
Fermi energy, a characteristic feature for the Fermi liquid. 

PACS numbers: 71.10.Fd, 71.27.+a 
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High T c superconductors found in cuprates have ex- 
hibited a rich phase diagram in doping (x = 1 — n) and 
temperature (T) plane [EJ. At half-filling the materials 
are antiferromagnetic insulators over a wide range of tem- 
perature up to 200-300 K, and away from half- filling an- 
tiferromagnetic long-range order is destroyed by 2-3 % 
doping and further doping makes the cuprates be metallic 
and eventually superconducting. Right after the discov- 
ery of high T c superconductors, the Hubbard model was 
proposed by Anderson 0] as the simplest model which 
might be able to describe the rich and anomalous phase 
diagram in the copper oxides. Although it was solved 
exactly in one-dimension the exact solution in more 
than 2-dimensions is not known yet. Nontherless, there is 
increasing numerical evidence j^H that for a half-filled 
2D band, the T = ground state is antiferromagnetic 
insulating for all U, while at finite temperatures strong 
fluctuations destroy long-range antiferromagnetic order 
due to the Mermin- Wagner theorem || . Away from half- 
filling, however, there have been lack of reliable studies on 
the systematic change of dynamical properties by doping 
which include, particularly, the low frequency behavior 
of the self-energy. 

Although quantum Monte Carlo (QMC) calculations 
have exhibited the systematic change of the spec- 
tral properties away from half-filling in the 2D Hubbard 
model, the lack of a reliable method for numerical ana- 
lytical continuation from QMC data containing inherent 
statistical errors as well as finite size effect have hindered 
a detailed study of some dynamical quantities such as 
the single particle self-energy. While there have been 
extensive studies on the low frequency behavior of the 
self-energy within fluctuation exchange (FLEX) approxi- 
mation @, 2D critical spin fluctuations are not properly 
taken into account in FLEX approximation. As a result, 
the scattering rates are linear (in frequency) over a wide 
range of frequency, but they always vanish quadratically 
near the Fermi energy for any available doping concen- 



trations within FLEX. 

Recently we formulated a theory to the 2D Hubbard 
model in a manner free of finite size effect and numerical 
analytical continuation, yet containing the essential fea- 
tures of the 2D Hubbard model, i.e., the correct atomic 
limit at large w and 2D spin fluctuations jioj . In this pre- 
vious work for a half-filled 2D band, 2D critical spin fluc- 
tuations are found to give rise to a strong local maximum 
at the Fermi energy in the scattering rates for low temper- 
atures, leading to a pseudogap in the spectral function. 
At zero temperature they are expected to yield eventu- 
ally an antiferromagnetic insulating state. Therefore, an 
important and natural question away from half-filling is 
whether strong 2D critical spin fluctuations are destroyed 
by doping or not, and if it is the case, how and at which 
doping concentration electrons lose their strong correla- 
tions, eventually leading to a Fermi liquid- like behavior. 
In this Letter, we would like to answer this one of the 
central issues in the high T c superconductivity for the 
first time. 

In order to take into account strong 2D critical fluc- 
tuations properly, we impose the following three exact 
sumrules to the spin, charge, and particle-particle sus- 
ceptibilities 
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^2xs P (q) =n- 2(n T nx) 

^2 Xch(q) =n + 2{n ] n i ) - n 2 
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T and N are the absolute temperature and number of 
lattice sites, q is a compact notation for (q, iu n ) where 
iv n are either Fermionic or Bosonic Matsubara frequen- 
cies. The dynamical spin, charge and particle-particle 
susceptibilities are calculated by 
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Xsp(q) = 

Xch(q) = 

Xpp(q) = 
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an d Xpp(?) are irreducible particle-hole and 
particle-particle susceptibilities, respectively, which are 
computed from 
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where G (k) is the noninteracting Green's function. 
U sp , U c h, and U pp in Eq. are renormalized inter- 
action constants for each channel which are calcu- 
lated self-consistently by making an ansatz U sp = 

^(«T«l)/((«T)(«i)) in E q- 0- B y defining U sp , U ch , 
and U pp this way, the Mermin- Wagner theorem as well as 
correct atomic limit for large uj are satisfied simultane- 
ously JTo| . In order to find the chemical potential for in- 
teracting electrons, first we calculate Eqs. (0)-(0) and the 
self-energy (Eqs. (0) in Ref. ]10[] ) with the noninteract- 
ing Green's function whose noninteracting chemical po- 
tential gives a desired electron concentration. Then, the 
chemical potential for interacting electrons is determined 
in such a way that the calculated electron concentration 
with the interacting Green's function is the same as the 
desired value. Throughout the calculations the unit of 
energy is t and all energies are measured from the chemi- 
cal potential /i. We used a 128 x 128 lattice in momentum 
space and performed the calculations by means of well- 
established fast Fourier transforms (FFT). It should be 
also noted that we used a real frequency formulation in 
Eqs. (0)-(0) to avoid any possible uncertainties associated 
with numerical analytical continuation. 

We start in Fig. [l] by studying the single particle spec- 
tral functions for U = 4 at n = 1.0, T = 0.05 (Fig. 0(a)) 
and at n = 0.87, T = 0.005 (Fig. 0(b)). Within our nu- 
merical calculations, we do not have the self-consistent 
solution for U sp at lower temperatures than 0.05 near 
half-filling. Note that the spin-spin correlation length 
exponentially increases like £ ~ exp(constant/T) at low 
temperatures near half-filling. In Fig. 0(a) for a half-filled 
band, an antiferromagnetic pseudogap is open at the 
Fermi energy along the antiferromagnetic zone boundary 
and a small shadow structure is found in the opposite 
frequency side of the main peak at the other momenta. 
At n = 0.87 (Fig. 0(b)), generally the spectral weight 
is transferred near the Fermi energy and a shadow fea- 
ture is strongly suppressed throughout the Brillouin zone. 
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FIG. 1. Spectral functions for U = 4 (a) at n = 1.0, 
T = 0.05 and (b) at n = 0.87, T = 0.005. 



The spectral functions along the antiferromagnetic zone 
boundary are so drastically changed at this electron con- 
centration that the antiferromagnetic pseudogap is com- 
pletely closed and a sharp quasiparticle peak begins to 
develop near the Fermi energy. 

In Fig. we present the systematic change of the den- 
sity of states and low frequency part of the self-energy 
by varying the electron concentration from n = 1.0 to 
0.95, 0.90 and 0.87. As the concentration is decreased, 
the spectral weight tends to move near the Fermi energy 
and to fill the antiferromagnetic pseudogap (Fig. 0(a)). 
At n = 0.90 (dotted curve), most of the antiferromag- 
netic pseudogap is closed and a small peak develops at 
the Fermi energy. At n = 0.87 (solid curve), the antifer- 
romagnetic pseudogap is completely destroyed by dop- 
ing with the appearance of a narrow and sharp spectral 
weight at the Fermi energy. At this stage a Fermi liquid- 
like behavior is expected to set in. In Fig. 0(b) the low 
frequency part of the scattering rates is shown for the 
same parameters as in Fig. 0(a). At n = 1.0 (dot-dashed 
curve), strong 2D critical spin fluctuations give rise to 
enormous scattering rates near the Fermi energy, even- 
tually leading to an antiferromagnetic insulator at zero 
temperature. At n = 0.95 (dashed curve), they are still 
persisting and shifted to the positive frequency side. At 
n = 0.90 (dotted curve), most of the critical fluctuations 
are destroyed by doping with a small remnant in the pos- 
itive frequency side. At n = 0.87 (solid curve), the 2D 
critical fluctuations are completely washed out and the 
imaginary part of the self-energy vanishes quadratically 
in frequency near the Fermi energy, a characteristic fea- 
ture for the Fermi liquid. 

In Fig. 0the single-particle spectral functions for U = 8 
at n = 1.0, T = 0.05 (Fig. |(a)) and at n = 0.80, 
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FIG. 2. (a) Density of states and (b) low frequency domain 
of the imaginary part of the self-energy at the noninteract- 
ing Fermi surface for U — 4 at different electron concentra- 
tions n = 1.0, 0.95, 0.90, and 0.87 denoted as the dot-dashed, 
dashed, dotted, and solid curves, respectively. For n = 1.0 
and 0.95, T is 0.05, while for n = 0.90 and 0.87, T is 0.005. 

T = 0.005 (Fig. |(b)) are shown. At n = 1.0, four peaks 
are found in the spectral function, two peaks associated 
with the antiferromagnetic bands near the Fermi energy 
and the other two with the Hubbard bands in the inter- 
mediate frequency regime. The Hubbard bands appear 
dispersionless since they arise from strong local repul- 
sions, and for U = 8 they are more pronounced than 
the antiferromagnetic bands JlO) ] . The antiferromagnetic 
pseudogap is more visible along the antiferromagnetic 
zone boundary, but its total spectral weight is much more 
suppressed compared with that for U = 4 in Fig. At 
n = 0.80 (Fig. ||(b)), the spectral weight is transferred 
near the Fermi energy throughout the Brillouin zone, as 
is the case for U = 4 in Fig. |](b) . At this electron concen- 
tration, the Hubbard bands are strongly suppressed and 
their significant spectral weight moves near the Fermi en- 
ergy to fill the antiferromagnetic pseudogap. The sharp 
quasiparticle peaks appear near the Fermi energy. 

The density of states and low frequency part of the 
scattering rates for U = 8 at different electron concen- 
trations from n = 1.0 to 0.91, 0.85 and 0.80 are presented 
in Fig. [|. As the concentration is decreased, a consider- 
able spectral weight associated with the Hubbard bands 
(Fig. |^(a)) moves near the Fermi energy and starts to 
fill the antiferromagnetic pseudogap. At n = 0.85 (dot- 
ted curve), most of the antiferromagnetic pseudogap is 
closed and a small peak develops at the Fermi energy, 
as in Fig. ||(a). At n = 0.80 (solid curve), the antiferro- 
magnetic pseudogap is completely closed by doping and 
a narrow, sharp spectral weight appears at the Fermi 
energy. In Fig. f|(b) the low frequency part of the scat- 



(0,0) 



J 




1 


(a) 


J 



































20.0 



0.0 



20.0 




-20.0 



20.0 



FIG. 3. Spectral functions for U = 8 (a) at n = 1.0, 
T = 0.05 and (b) at n = 0.80, T = 0.005. 



tering rates is also shown for the same parameters as in 
Fig. Q(a). At n — 1.0 and 0.91 (dot-dashed and dotted 
curves), strong 2D critical fluctuations are still dominat- 
ing near the Fermi energy. At n — 0.85 (dotted curve), 
most of the 2D critical fluctuations are already destroyed 
by doping with a small remnant. Eventually 20 % doping 
(solid curve) destroys the 2D critical fluctuations com- 
pletely and the scattering rates approach zero quadrati- 
cally in frequency near the Fermi energy, as in Fig. ||(b) 
for (7 = 4. 

We plot in Fig. |5| the imaginary parts of the dynami- 
cal spin susceptibility for U = 8 at n = 1.0, T = 0.05 
(Fig. |(a)) and at n = 0.80, T = 0.005 (Fig. |(b)). 
Both cases in Fig. || show a similar behavior at small 
and intermediate momenta except that the maxima for 
the former occur at a smaller energy scale than that for 
the latter, consistent with the recent QMC calculations 
S. The most dramatic difference happens at q = (7r, tt) 
(solid curves). For a half-filled 2D band (See the inset in 
Fig. ||(a)), Imxsp{q-,v) shows a nearly divergent behav- 
ior at small frequencies (its maximum 6760/i occurs at 
5.5 x 10~ 5 i), while at n = 0.80 (Fig. |(b)) Im Xsv {q,v) 
exhibits only a broad maximum at 0.75<. This drastic 
change in Imxsp (<f, v) is mainly responsible for the ma- 
jor differences discussed earlier for these two electron con- 
centrations. 

Before closing we comment on the low frequency be- 
havior of the self-energy. Iml](k, oj) for U = 4, n = 0.87 
and for U = 8, n = 0.80 are also examined on a much 
smaller energy scale from -0.1 to 0.1 which is not shown 
in this Letter. The quadratic behavior is still found with 
a small constant shift (due to a finite temperature effect). 
For the former case, the scattering rates vanish as ui 1 - 95 
near the Fermi energy, while for the latter as w 1 ' 94 . As 
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FIG. 4. (a) Density of states and (b) low frequency domain 
of the imaginary part of the self-energy at the noninteract- 
ing Fermi surface for U — 8 at different electron concentra- 
tions n = 1.0, 0.91, 0.85, and 0.80 denoted as the dot-dashed, 
dashed, dotted, and solid curves, respectively. For n = 1.0 
and 0.91, T is 0.05, while for n = 0.85 and 0.80, T is 0.005. 

is shown in Figs. ^ and |], we do not also find any indica- 
tion of a linear frequency dependence for the scattering 
rates. Hence, as long as the 2D Hubbard model with only 
nearest neighbor hopping is concerned, our calculations 
do not show the frequency dependence of the scattering 
rates from FLEX approximation Q or from the marginal 
Fermi liquid hypothesis 

In summary, the doping dependence of the spectral 
functions, density of states and low frequency behav- 
ior of the self-energy has been studied for the 2D Hub- 
bard model based on our recently developed theory for 
the Hubbard model. Strong 2D critical spin fluctuations 
dominating near half-filling are completely destroyed by 
13 and 20 % doping concentrations for U = 4 and 8, 
respectively. Below these concentrations the scattering 
rates vanish quadratically in frequency near the Fermi 
energy, a characteristic feature for the Fermi liquid. 
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FIG. 5. Imaginary part of the dynamical spin susceptibility 
for U = 8 (a) at n = 1.0, T = 0.05 and (b) at n = 0.80, 
T = 0.005. The dot-dashed, dashed, dotted and solid curves 
denote ImXs P (q,v) at i = 1, 2, 3 and 4 for q — i(ir/A, 7r/4), 
respectively. 
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